Phase transition and critical behavior of d=3 chiral fermion models with left/right 

asymmetry 
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We investigate the critical behavior of three-dimensional relativistic fermion models with a 
U(-^l)l X U(1)r chiral symmetry reminiscent of the Higgs- Yukawa sector of the standard model 
of particle physics. We classify all possible four-fermion interaction terms and the corresponding 
discrete symmetries. For sufficiently strong correlations in a scalar parity-conserving channel, the 
system can undergo a second-order phase transition to a chiral-symmetry broken phase which is 
a 3d analog of the electroweak phase transition. We determine the critical behavior of this phase 
transition in terms of the critical exponent v and the fermion and scalar anomalous dimensions for 
A^L > 1. Our models define new universality classes that can serve as prototypes for studies of 
strongly correlated chiral fermions. 



I. INTRODUCTION 

The interplay between statistical physics and particle 
physics has stimulated substantial progress in both areas 
in the last decades. It has given rise to the renormaliza- 
tion group (RG) [I!, S i which has led to a deep under- 
standing of critical phenomena on the statistical-physics 
side and of the possible structure of fundamental inter- 
actions on the particle-physics side. In most cases, this 
interplay arises from analogies on the technical level of 
dealing with fluctuations either of statistical or quantum 
nature. 

In the present work, we intend to push this interplay 
one step further: we aim at a construction of a statistical- 
physics model (or a field-theoretic variant thereof) which 
has structural similarities to a crucial building block of 
the standard model of particle physics: the Higgs- Yukawa 
sector. We are motivated by the fact that a profound 
understanding of electroweak symmetry breaking in the 
standard model or alternative scenarios may require a 
quantitative control of fluctuating chiral fermions and 
bosons beyond perturbation theory.^ This is a challeng- 
ing problem which is even lacking a simpler benchmark 
example where nonperturbative techniques can prove 
their validity. 

A stage for quantitative comparisons between field- 
theoretical tools is set by critical phenomena of lower- 
dimensional systems, most prominently the computation 
of critical exponents of second-order phase transitions. 
In fact, relativistic fermionic models such as the Gross- 
Neveu model in 3d are known to exhibit a second-order 
phase transition; corresponding studies of the critical be- 
havior have been performed by various methods 0]- 



^ This work has partly been inspired by a recently discovered 
asymptotic-safety scenario for the Higgs- Yukawa sector of the 
standard model which aims at a solution of the triviality and 
an improvement of the hierarchy problem of the standard-model 
Higgs sector 



Near the phase transition, the critical behavior is de- 
termined by fluctuations of fermions as well as bosonic 
bi-fermion composites; the expectation value of the latter 
also serves as an order parameter for the phase transition. 

In the same spirit, we devote this work to a construc- 
tion of chiral four-fermion models with a U(A^l)l x U(1)r, 
symmetry. Such models support a left/right asymmetry 
similar to the chiral structure of the Higgs- Yukawa sector. 
If such models undergo a phase transition to an ordered 
phase with broken chiral symmetry this transition can be 
viewed as an analogue of the electroweak phase transition 
in 3d. The corresponding critical behavior defines new 
universality classes, the properties of which are a pure 
prediction of the system. 

In the past, d = 3 dimensional fermionic systems with 
left /right symmetric chiral symmetries such as QED3 or 
the Thirring model have been under investigation in a va- 
riety of scenarios [1, S [13, Ell, [13, with apphcations to 
condensed-matter physics, high- Tc cuprate superconduc- 
tors and, recently, graphene [31 . In some of these 
models, the number of fermion flavors serves as a control 
parameter for a quantum phase transition. As the critical 
number of fermions is an important quantity, nonpertur- 
bative information about these models for varying flavor 
number A^f is required. Since good chiral fermion proper- 
ties for arbitrary TVf for massless fermions still represents 
a challenge, e.g., for lattice simulations, other powerful 
nonperturbative techniques are urgently needed. 

In this work, we investigate this class of chirally sym- 
metric fermion models with left /right asymmetry by 
means of the functional RG. The functional RG and as- 
sociated RG flow equations such as the Wetterich equa- 
tion [l6l | are an appropriate tool for this purpose, since 
fermions and bosons can both be treated equally well. 
Also, the flow equation is an exact equation and facili- 
tates the construction of non-perturbative approximation 
schemes, see for reviews. 

After a classiflcation of possible fermionic models in 
Sec. HTJ we concentrate on strong correlations in a scalar 
parity-conserving channel. We study the critical behav- 
ior of a possible condensation of bosonic bi-fermion com- 
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posites into this channel. The resulting effective scalar- 
fermion Yukawa model is reminiscent to the standard- 
model Higgs- Yukawa sector. We analyze this effective 
Yukawa model in Sec. IIIII with the aid of the functional 
RG, deriving the nonperturbative RG flow equations to 
next-to-leading order in a derivative expansion. 

The left /right asymmetry of our model is controlled by 
varying the number A^l of left-handed fermions. This im- 
balance provides us with an external control parameter 
for the relative amplitude of boson and fermion fluctua- 
tions. This allows us to vary the fixed-point properties 
of the RG flow, implying also a variation of the critical 
properties of these systems. In particular, the fixed-point 
potential for the bosonic order parameter can be in the 
symmetric regime for small Ni^ or in the broken regime 
for larger A^l, as discussed in Sec. IIVI The small A^l 
regime turns out to be particularly interesting, as the 
anomalous dimensions of all fields have to satisfy a sum 
rule within our truncation, in order to give rise to a non- 
trivial fixed point in the Yukawa coupling. 

The quantitative reliability of our results for the crit- 
ical properties can be checked by a reduction of our 
model to a corresponding purely bosonic O(A^) model. 
For the latter, our method has exhaustively been in- 
vestigated and used with a remarkable quantitative suc- 
cess [ITI . [isl . IT9| : the results for critical exponents have 
reached an accuracy which is comparable if not superior 
to that of, e. g., l attice simulations and high-order e ex- 
pansions, see [l^ for a recent comparison. Our inclusion 
of fermions can reliably be based on this footing, as our 
models approach an 0(2iVL) model in the limit of large 
TVl; this is because the composite scalar degrees of free- 
dom dominate the fluctuation contributions as a result 
of the chiral structure. For smaller A^l, our results show 
quantitative as well as qualitative differences to 0{N) 
models, as expected due to the fermionic contributions. 
The resulting models and quantitative findings consti- 
tute and characterize a new set of universality classes, 
classified by the chiral symmetry content. In particular, 
we provide for quantitative predictions for the critical 
exponents for these universality classes, which have not 
been investigated with any other method so far. We be- 
lieve that these can serve as a first benchmark for other 
nonperturbative methods which are urgently needed for 
a study of chiral phase transitions in strongly correlated 
chiral fermions. 



II. CLASSICAL ACTION AND SYMMETRY 
TRANSFORMATIONS 

Let us first consider a general fermionic model in 
d = 2 + I Euclidean dimensions with local quartic 
self- interaction, being invariant under chiral U(A^l)l ® 
V{^r)r transformations. The Dirac algebra 

{7m' } = 2(5^1^, (1) 



could minimally be realized by an irreducible represen- 
tation in terms of 2 X 2 matrices. As this representation 
does not permit a chiral symmetry, the massless theory 
could not be separated from the massive theory by an 
order-disorder transition. We therefore work exclusively 
with a 4 X 4 reducible representation of the Dirac algebra 

^^=(i', >)' ^'=^'2,3, (2) 

with {(Ti, a2, era} being the 2x2 Pauli matrices. Such 
models with four-component fermions in three dimen- 
sions have extensively been studied in the past (e.g. in 
the context of spontaneous chiral symmetry breaking in 
QED3 or the Thirring model) There has been 

renewed interest in the last years, due to potential appli- 
cations to high-Tc superconductivity in cuprates or elec- 
tronic properties of graphene [HI, [13, ■ For a review, 
see e.g., |13| . There are now two other 4x4 matrices 
which anticommute with all 7^ as well as with each other, 

(1 0) ^^'^ ^ 71727374 = (^Q ■ (3) 

Together with 

1, CT^i^ — ^[7^,7./] ifJ- < 17^.74, 17^75, 17475, (4) 

these 16 matrices form a complete basis of the 4x4 Dirac 
algebra, 

{7a}a=i 16 = { 1,7m, 74, crM'',i7/x 74, 17^75, 17475, 75 } ■ 

(5) 

Now, we define chiral projectors 

PL/R = i(l±75) (6) 

that allow us to decompose a Dirac fermion i/i into the 
left- and right-handed Weyl spinors ipi^ and 'i/'R, 

V'L/R = P\./Ki\ V'L/R = V'pR/L- (7) 

('0 and Tp are considered as independent field variables 
in our Euclidean formulation.) Note that there is a cer- 
tain freedom of choice of the notion of chirality here: 
we could have chosen just as well Pl/r = (1 ± 74)/2 or 

^L/R = (1 ±i7475)/2 as chiral projectors. This would 
have led us to different definitions of the decomposi- 
tion into Weyl spinors. All these chiralities remain con- 
served under Lorentz transformations since all three pro- 
jectors commute with the generators of the Lorentz trans- 
formation of the Dirac spinors, [75,cr^i/] = [74,cr/ji/] = 
[17475,^^1.] = 0. 

We consider N-r, right-handed and A^l left-handed 
fermions, where N^i and A^l do not have to be identi- 
cal. We impose a chiral 'U{N\,)\, ® U(A^r)r symmetry 
with corresponding field transformations which act inde- 
pendently on left- and right-handed spinors. 

U(ArL)L: 0£-C/£V^ ri.^i'iiUlr. (8) 
U(A^r)r: ^R-^^K)''^ (9) 
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Here, C/l and JJr, are unitary A'l x TVl and A^r x TVr ma- 
trices, respectively. For = e'" 6°-'' and C/^^ = e"'" S'''' 
we obtain the usual U(1)a axial transformations, whereas 
for Ul'' = e'" d"'' and U^^ = e'" 6"'' we get U(l)v phase 
rotations. The symmetry thus is 

U(7VL)L®U(iVR)R, = SU(ArL)L®SU(7VR)R®U(l)A®U(l)v, 

(10) 

with chiral SU(AfL.R) factors. 

Due to the reducible representation of the Dirac alge- 
bra, there is also some freedom in the definition of the 
discrete transformations [lo| . Charge conjugation may 
be implemented by cither 

C ■■ ^l/K ^ {^I/rCY - ^t/R - - (C^-'WR,)^ - (11) 

with C ~ 7275, or 

C ■■ r^/R - {rn/Lcf , ^L/R - - , (12) 

with C ~ 7274, or a unitary combination thereof. In the 
same manner, the parity transformation corresponding 
to 

{xi,X2,X3) 1-^ i-Xi,X2,X3) =: X, (13) 

with {xi,X2) as space coordinates and x^ as the (Eu- 
clidean) time coordinate, may be implemented by either 

^ : i'l/Ri^) ^ ^^Wr(5), i'l/Ri^) ^ i'l/Rii)P^, 

(14) 

with P ~ 7i74, or 

'P ■■ ^'IM^) ^ P^R/Lii), ^L/r(^) ^ i'R/Lii)P\ 

. ^^^^ 

with P ~ 7i75. Similarly, time reversal corresponding to 

{xi,X2,X3) 1-^ {xi,X2,-X3) X (16) 

reads cither 

r : ^l/^ix) ^ T^£/r(x), V5£/j,(x) ^ V;L/R(i)rt, (17) 
with T ~ 7273, or 
r : ^-/i,(x) f^g./L(x), ^Pl^^{x) ^ V;S,/L(i)f t, (18) 

with f = 71. Note that every matrix e {C, C, P, P, T, f} 
is unitary, such that charge conjugation and parity inver- 
sion are unitary, and time reversal is anti-unitary. 



In order to derive the explicit transformation proper- 
ties of the bilinears, it is useful to recall that 71 and 73 
are antisymmetric and purely imaginary, whereas 72, 74, 
and 75 are symmetric and real. The results are listed in 



Table |T1 where we have introduced 

7m — (-71,72,73)/., (19) 

(o'l2, CT13, 0'23) (-0-12, — (T13, (T23), (20) 

7m := (71,72,-73)/., (21) 

(o'i2, CT13, o'23) := (cti2, — CT13, — (T23). (22) 



These transformation properties facilitate a discussion of 
possible bilinears and 4-fermi terms in the action of our 
model. In addition to Lorentz invariance, we impose an 
invariance of our theory under U(A^l)l ® U(A^r)r chiral 
transformations, C charge conjugation, V parity inver- 
sion, and T time reversal. The theory then automatically 
is also invariant under CV, VT , and CT transformations, 
since CV ^CV,Vf ^ VT, and Cf = CT. (The equiva- 
lence holds up to U(l)v phase rotations of the spinors.) 
As a consequence, no bilinears to zeroth order in deriva- 
tives are permitted. To first order, only the standard 
chiral kinetic terms 

i^lid^-f^^l and V;Ri5/.7/.V'R (23) 

can appear. In particular, all possible mass terms are 
excluded by symmetry: V'lV'r and V'rV'l ^'''^ '^ot chirally 
symmetric, and '0l74''/'l ^gW. as '0r74V'r are not invari- 
ant under V parity transformations. The same holds for 
terms involving 17475. On the level of 4-fermi operators, 
the interaction terms must have the form 

{4^iiAi^\) {^iiAi^D , {rRiAi^i) (vi^7AV'R) (24) 

with 7A e {7m>74}, 

i^llB^li) {i^RlB^D with 7Be{l,i7/.74}, (25) 
or, with inverse flavor structure, 

{^rIa^r) (^lTaV'l) • (26) 

Terms with 7^ S {i7/. 75, 17475} or 75 e {75,0-^1/} are 
equal to these up to a possible sign, since ■(/'l and i/jr are 
eigenvectors of 75, and cr^i, = —'^c^vplplAlb- Terms with 
lA e {1,75, 17/^74, cr^!/} or 7B e {7/x, 17^75, 74, 17475} are 
identically zero, since P^Pl = PlPr = 0. The terms in 
Eq. are not independent of the terms in Eqs. (|24p . 
P5|l . but are related by Ficrz transformations: 

i^h^i^l) + 1 (^£74^^) (V5^4V'£) , (27) 

i^h^^^'t) - \ (^£740 (v5^4V'£) , (28) 
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TABLE I: Properties of fermion bilinears under discrete transformations. The arguments of the transformed fields are x = 
{—Xi,X2,x-i) in the case of parity and x — {xi, X2, —xs) in the case of time reversal. The bilinears with (L ^ R) transform 
analogously. 
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V'lc^m.-V'r 
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-^r74V'r 


-^£74V'£ 




^r74^^ 


V'£74V'£ 


-^r74V'r 


V'£i7/.74^/'R 


V'Ri7K74^£ 


-^£i7M74i/'R 


-V'£i7M74^R 




^Ri7M74V'£ 


V'£i7M74V'R 


-V'Ri7M74< 








^ (V^rV^l) (^lV-r) - 


1 
2 


(^Ri7M74^^) 


(V;^7m74V'r) , 


(29) 






ii'llii'l) = - 


^ (V^rV^D + 


1 
2 


(^Ri7M74^^) 


(V;^7m74V'r) . 


(30) 



Here, we have suppressed the analogous equations with 
(L <-!■ R) for simplicity. We thus end up with six indepen- 
dent 4-fermi terms preserving U(A^l)l ® ^{N-r)k chiral 
and C,V, and T symmetry, 

(^£^r) {i'ii^D , (31) 

(^^£74^^) {i^ilAi^D , {^Mi) (474^r) , (32) 

(^£i7M74V'R) (^Ri7M74^£) . (34) 

Note that the corresponding set in d = 3 + 1 dimensions 
would be smaller as the 74 terms would be constrained 
by a larger Lorentz symmetry. 

In a (partially) bosonized language after a Hubbard- 
Stratonovich transformation, we encounter six boson- 
fermion interactions: The first one corresponding to 
Eq. ((3T|) couples the fermions to a scalar boson (scalar 
with respect to V parity), the second and the third ([5^ 
to a pseudo-scalar boson, the fourth and the fifth ([55]) 
to a vector boson, and the sixth ([M)) to a pseudo-vector 
boson. Further bosonic structures in the flavor-singlet 
channels appear in the corresponding Ficrz transforms 
of Eqs. 

We expect that a general model based on these interac- 
tions exhibits a rich phase structure, being controlled by 
the relative strength of the various interaction channels. 
Aiming at an analogue of the electroweak phase transi- 
tion, we focus in this work on the V parity-conserving 
and Lorcntz-invariant condensation channel, parameter- 
ized in terms of the first interaction term. This channel 
is also invariant under T time reversal, whereas the bo- 
son transforms into its complex conjugate under C charge 
conjugation. Moreover, we confine ourselves to the case 
of A'r = 1 right-handed fermion flavor and A^l > 1 left- 
handed fermion flavors, allowing for a left /right asym- 
metry similar to the standard model of particle physics 
(where A'l = 2). The microscopic action of our model in 



the purely fermionic language then reads 

54-fcrmi = y"d=^a;{Vi£i^V'£ + ^Ri^V'R 

+2A(^£^r)(^r^£)}. (35) 

Via Hubbard-Stratonovich transformation, we obtain the 
equivalent Yukawa action 

+0"V"r^£-0V£V'r}, (36) 

where the complex scalar serves as an auxiliary field. 
The purely fermionic model can be recovered by use of 
the algebraic equations of motion for and (t>'^''\ 

0"=-2AV5rV£, ^^=2X4,1^;^. (37) 

Here, we can read off the transformation properties of 
the scalar field under the chiral symmetry, 

^« ^ Ul''(j)''Ul, (t)"'' ^ Uncl)''Hul)'"' ■ (38) 

The composite scalar field 0" represents an order param- 
eter for an order-disorder transition. As long as 0" has 
a vanishing expectation value, the system is in the sym- 
metric phase with full chiral V{Ni^)-l x U(1)r symmetry; 
the fermions are massless, whereas the scalars are generi- 
cally massive as determined by the symmetry-preserving 
effective potential for the scalars. If 0" acquires a vac- 
uum expectation value the chiral SU(AfL) factor is broken 
down to a residual SU(AfL — 1) symmetry. In addition, 
the axial U(1)a is broken, whereas the charge-conserving 
vector U(1)v is preserved. In the broken phase, the spec- 
trum consists of one massive Dirac fermion, one massive 
bosonic radial mode, A^l — 1 massless left-handed Weyl 
fermions and 2A^l — 1 massless Goldstone bosons. 

Near the phase transition, we expect the order- 
parameter fluctuations to dominate the critical behavior 
of the system. Universality suggests that the degrees of 
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freedom parameterized by the action (j36p are sufficient to 
quantify the critical behavior of this transition, indepen- 
dently of the presence of further microscopic fermionic 
interactions of Eqs. (PT|) - ([M)) . Concentrating on the ac- 
tion ([36]), we observe that the purely scalar sector, i.e., 
the scalar mass term, has a larger symmetry group of 
0(2A'L)-type. It is therefore instructive to compare the 
critical behavior of our fermionic model with that of a 
standard scalar 0{2Ni^) model which is known to un- 
dergo a second order phase transition associated with 
a Wilson-Fisher fixed point. Differences in the corre- 
sponding critical behaviors can then fully be attributed 
to fermionic fluctuations near the phase transition. 



III. EFFECTIVE AVERAGE ACTION AND RG 
FLOW 

Integrating out fluctuations momentum shell by mo- 
mentum shell near the phase transition described above, 
the effective action (effective Wilson-type Hamiltonian) is 
expected to acquire all possible operators of mixed scalar 
and fermionic nature compatible with the symmetries of 
the action (j36p . In the present work, we constrain the 
flow of this action functional to lie in the subspace of the 
full theory space spanned by the following ansatz valid 
at an RG scale k: 



Tk = 



h^iri^n}- (39) 



The fermion fields tpf^ and ipR have standard kinetic terms 
but can pick up different wave function renormalizations 
Zh.k and ZR.fc. The index a runs from 1 to A^l- The 
bosonic sector involves a standard kinetic term with wave 
function renormalization Z^^k and an effective potential 
Uk{p), where p = It is sometimes useful, to ex- 

press the complex scalar field in terms of a real field basis 
by defining 

,^« = -i=(0?+i<^^), 0-t = _L(0j_i0^-), (40) 

where (/)i,(/'2 € R. AH parameters in the effective aver- 
age action are understood to be scale dependent, which 
is indicated by the momentum-scale index k. The scale 
dependence is governed by the Wetterich equation, which 
allows for a nonperturbative construction of quantum 
field theory in terms of the effective average action Tk 



9,rfc[$] = iSTr^ + 



(dtRk) 



(41) 



(2) 

Here, Tj. ' [$] is the second functional derivative with re- 
spect to the field <!>, the latter representing a collective 



field variable for all bosonic or fermionic degrees of free- 
dom, and Rk denotes a momentum-dependent regulator 
function that suppresses IR modes below a momentum 
scale k. The solution to the Wetterich equation provides 
an RG trajectory in theory space, interpolating between 
the bare action iSvuk to be quantized Ffe^A S'vuk and 
the full quantum effective action F = Ffc^o, being the 
generating functional of IPI correlation functions. For 
reviews, see e.g., [l7j . 

The ansatz ([5^ . in fact, represents the next-to- leading 
order in a systematic derivative expansion of the effective 
potential in the scalar sector and a leading-order vertex 
expansion in the fermionic sector. (The leading-order 
derivative expansion is obtained by setting all wave func- 
tion renormalizations 2^(</).L.R),fc = const.) This expansion 
can consistently be extended to higher orders and thus 
defines a legitimate and controllable nonperturbative ap- 
proximation scheme. It has indeed proved its quantita- 
tive reliability already in a number of examples involving 
Yukawa sectors 0, Ho, [H HI . 

In order to fix the standard RG invariance of field 
rescalings, we define the renormalized fields as 



= Z 



1/2 J, 7 '7I/2 ; 



(42) 



For the search for a fixed point where the system is 
scale invariant, it is useful to introduce dimensionless 
renormalized quantities. In order to display the dimen- 
sion dependence, we perform the analysis in d spacetime 
dimensions, where the dimensionless renormalized field, 
Yukawa coupling and scalar potential read 



Uk{p) = k Ukip)\p=k''-'2p/z^,,- 



(43) 
(44) 
(45) 



The flow of the wave function renormalizations Z^,k, ■^L,fc 
and Zfi fc can be expressed in terms of scale-dependent 
anomalous dimensions 

770 = -dt\nZ^,k, ?7l/r = -9tlnZL/R„fc . (46) 

With these preliminaries, we can compute the flow of the 
effective potential (see App. IX)) . 

dtuk = -duk + p{d-2 + ri^)u'k + 2vd {2Nl - l)/o("fc) 



(47) 



where ^ = 2'*+i7r'*/2r(d/2), d.^ = 4 is the dimension of 
the 7 matrices, and the threshold functions. 



,(F)rf 
n,L/R 



2((5„,o + "-) 



2((5„,o + n) 



1 - 



1 



V4> 

d + 2) (l + cj)"+i' 
?7l/r \ 1 



\ d+lJ 



d+lJ (1 + cj)"+ 



T^48) 
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encode the information about a possible decoupling of 
massive modes. For the momentum regularization, we 
have used a linear regulator function j?^ here which is 
optimized for the derivative truncation [23| . 

Whereas the flow of the Yukawa coupling is unam- 
biguous in the symmetric regime, the Goldstone and 
radial modes can generally develop different couplings 
in the broken regime. Here, we concentrate on the 
Goldstone-mode Yukawa coupling to the fermions, as the 
radial mode becomes massive and decouples in the broken 
regime. In both phases, the flow of the Yukawa coupling 



can be written as (see App. [X]) 

dthl = iv^ + VL + VR + d - 'i)hl (49) 
-8vdhlKku'llf^f^'^{Kkhl,u'f^ + 2Kku'l,u'i,), 

where the scalar-potential terms have to be evaluated on 
the fc-dependent minimum of the potential pmin = i^k- 
In the symmetric regime, we have, of course, Hk ~ 0. 
The threshold function occurring in Eq. (|49p reads for 
the linear regulator: 



1 



d{l + cji)"i (1 + ws)"^ (1 + ws)"^ 



ni 



I +Ll>1 



d + l 



n-2 



1+^2 



'?0 

d+2 



"-3 



1 + CJ3 



V4> 
d + 2 



Setting the anomalous dimensions to zero defines the leading-order derivative expansion. At next-to- leading order, it 
is important to distinguish between Zi^^k and ZB.,k as they acquire different loop contributions, see below. The flows 
of the anomalous dimensions read (see App. [A| 



k"'-4 y'^kiik) 



8Wd,2r (FB)d.,2 / , r, ii\ : (FB)d,,2 t \\ 



(50) 



d 



"H2 [hk'ik:Ui, + 2KkUi,)+m\2 [hkKk,Uk) +2(Nl- 1)77112 {0,Uk) 

I 



where the potential is again evaluated at p = k^. Here, 
we have also introduced the corresponding threshold 
functions for the linear regulator 

I 

(F)d. X 1 

Too (UJ) = —r, 

TO4 ' (tj) 



(1+CJ)4 ' d-2 (l+Uj)^ 
(I - 77^ 1 \ 1 



\2d-A 4/ (1 +cj)2' 



^V-dT-l) (1+C.0-(1+C.2)-' 

(51) 

and rj^, := i(r/R, + t^l)- 



IV. FIXED POINTS AND CRITICAL 
EXPONENTS 

The Wetterich equation provides us with the flow of the 
generalized couplings dtgi = (3i{gi,g2, . . . ) of our trunca- 



tion, where the gi correspond to the Yukawa coupling 
and, e.g., the expansion coefficients of the scalar poten- 
tial, etc. A fixed point g* is defined by 



V« : Pt{gl,g2,.:) =0- 



(52) 



If the fixed point separates the disordered from an or- 
dered phase, it is a candidate for a second-order phase 
transition. For the investigation of the critical properties 
of the theory near this transition, we consider the fixed- 
point regime, where the flow can be linearized around the 
flxed point. 



dtgt ^ B^' {g* - gj) + . 



m 



a=g* 



(53) 



Let us denote the eigenvalues of the stability matrix Bf 
by w/. The index / labels the order of the eigenvalues 
according to their real part, starting with the smallest 
one which we call luq- Negative cui correspond to RG 
relevant directions for which an IR (UV) flxed point is 
repulsive (attractive). In analogy to the notion of critical 
phenomena in 0(A'^)-models, we deflne v = — I/wq, char- 
acterizing the critical exponent of the correlation length 
near the critical temperature. This is indeed justiflcd, 
as the largest critical exponent is associated with the 
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strongest RG relevant direction, being in turn related to 
the distance from the critical temperature. Thus, v, in 
fact, corresponds to the standard correlation- length ex- 
ponent. The subleading exponent is traditionally called 

w = Wi. 

In order to analyze the fixed point structure of this 
model, we will distinguish two different regimes of the 
system namely the symmetric regime (SYM), where the 
vacuum expectation value (vev) of the boson field is zero, 
and the regime of spontaneously broken symmetry (SSB), 
where it is nonzero. As boson and fermion fluctuations 
generically contribute with opposite sign, the existence of 
a fixed point requires a balancing between both contribu- 
tions together with potential dimensional scaling terms 
(such as, e.g., the first term on the right-hand side of 
Eq. We observe that this balancing is indeed 

possible in both regimes, depending on the number of 
left-handed fermion flavors. The origin of this flavor- 
number dependence is illustrated for the running of the 
scalar mass term or vacuum expectation value in Fig. [T] 
Whereas the scalar loop carries a weight ~ TVl as all 



pansion of the effective potential u 



A2 



FIG. 1: Loop contributions to the renormalization flow of 
the mass or the vev. The left loop involves only inner boson 
lines. The vertex A2 allows for a coupling between all available 
boson components. This implies a linear dependence on A^l 
for the renormalization of the boson contribution, see below. 
On the right panel, we depict the corresponding fermion loop 
contribution. The incoming boson 0a fully determines the 
structure of the fermion loop and does not allow for other left- 
handed inner fermions than ■(/)£, inhibiting an A^l dependence 
of the algebraic weight of this loop. 

scalar degrees of freedom can contribute, the left/right 
asymmetry structure leads to a weight ^ for the 

fermion loops in Fig. [TJ In the following, we analyze the 
different regimes in detail. 



A. The symmetric regime 

Here and in the following, we drop the subscript k indi- 
cating the scale dependence of the running couplings for 
simplicity. The presence of this scale dependence will be 
implicitly understood. Also we restrict our investigations 
to d = 3 from now on. 

In the symmetric regime, we employ the following ex- 



-TP 



(?7i^ = Ai). The minimum is assumed to be at p = 0, 
implying k = 0. This also implies that m? > 0. This 
expansion is plugged into the flow equations for the ef- 
fective potential, the Yukawa coupling and the anomalous 
dimensions. An expansion of the whole flow equation in 
terms of p yields the flows of the different running cou- 
plings A„. 

First, we write down the explicit expression for the 
Yukawa coupling, which in the SYM regime reduces to 



dth^ 



+ + m. - 



(55) 



This tells us that an interacting fixed point {h 7^ 0) can 
only occur if 



(56) 



Conversely, if a fixed point exists in the SYM regime, 
the sum rule ([55]) has to be satisfied by the anomalous 
dimensions. This statement holds exactly in the present 
truncation but may receive corrections from higher or- 
ders. We comment further on the relevance of this sum 
rule in the conclusions. The expressions for the anoma- 
lous dimensions read 



V<f> 



37r2 



(5 - ?/L 



vr), 
1 



(1 



1 



(57) 
(58) 

(59) 



This is a linear system of equations which can be solved 
analytically. Its solution expresses the anomalous dimen- 
sions in terms of the couplings and m^, 



2h^{2h^{Ni^ + 1) - 157r2(l + m^Y) 
h'^{Ni^ + 1) - 187r4(l + m2)2 

;i2(5/i2_ 127r2) 



1) - 187r4(H-m2)2^ 

/l2(5/i2_ 127^2) 



/i4(iVL -f- 1) - 187r4(l 



(60) 
(61) 
(62) 



These expressions can be plugged into the sum rule ([5l 
resulting in a conditional fixed point for depending on 



the size of 1 



which reads 



J 



/)2 

"■cond 



3^2 



8(iVL 



7 + 5m{2 



+ 2iVL - V33 + m(2 + m)(54 + 25m(2 + m)) + 12iVL -|- 4m(2 + m)Ni^ + ANl 
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For another solution with a positive root, we have not 
been able to identify a true fixed point of the full system 
by numerical means. Hence, this solution is ignored in 
the following. The solution with the negative root, how- 
ever, does give a fixed point and will be analyzed in the 
following. This solution is positive for all > and 
monotonously increasing. It ranges between h^{rn^ = 

0) = ^ (7 + 2iVL - y^33Tmj3TN^) /(^L + 1) 

and h?{m? = 00) = which is very convenient be- 
cause it ensures that the fixed point value for h? is 
bounded from above and from below in a very nar- 
row window for all . A true fixed point of the sys- 
tem requires fixed points for all scalar couphngs (m^(= 
Ai), A2, A3, ...). The flow equation of a coupling A,i is al- 
ways a function of the lower order couplings from the 
effective potential up to A„+i, i.e. 



C^tAn — fnih^, Ai, A„+i). 



(63) 



Inserting the conditional fixed point for the Yukawa cou- 
plings ([63|l into these flows leaves us with the problem of 
searching for a scalar fixed-point potential. This problem 
is familiar from scalar theories where the Wilson-Fisher 
fixed point follows from an equation similar to Eq. ([63]) 
(with = 0, of course). We solve the fixed-point equa- 
tions dtXn — approximately by a polynomial expan- 
sion of the potential up to some finite order n < rimax- 
Dropping the higher-order couplings Xn>n^^^^ ~ 0, the 
resulting system of fixed-point equations can be solved 
explicitly. We find suitable fixed-point solutions in the 
symmetric regime for Nj^ g {1,2}. The non- universal 
fixed-point values as well as the universal values for the 
anomalous dimensions at the fixed point and the first 
two critical exponents can be read off from Table HTl For 
these results, we have expanded the effective potential 
up to Ag at next-to-leading order in the derivative ex- 
pansion and computed the corresponding stability ma- 
trix, c.f. Eq. (|53p . including the Yukawa coupling flow. 
We would like to stress that the fixed-point equations in 
the present case are technically much more involved in 
comparison with those of scalar O(A^) models in a sim- 
ilar approximation, as the insertion of the conditional 
Yukawa coupling fixed point introduces a much higher 
degree of nonlinearity. 



nature of the phase transition and the corresponding crit- 
ical behavior is characteristic for our fermionic model. In 
particular for small TVl, the fermionic fluctuations con- 
tribute with a comparatively large weight to the critical 
behavior, as discussed in Fig. [T] 

In order to estimate the error on our results arising 
from the polynomial expansion of the effective potential, 
we study the convergence of the fixed-point values and of 
the critical exponents as a function of increasing trunca- 
tion order for TVl = 2. Figure [5] displays our results for a 
truncation beyond order p" ~ 0^" for n = 2, 4, 6, 7. All 
quantities show a satisfactory convergence with a varia- 
tion on the 1% level among the highest-order truncations. 
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FIG. 2: Error estimate for A^l = 2: fixed-point values 
, , A2 and critical exponents UjCj as a function of the 
highest-order term n in the p expansion of the effective poten- 
tial. Left panel: hi (blue dots), A2 (purple squares). Right 
panel: m* (red dots), v (green squares), lo (orange diamonds). 



TABLE XL Fixed-point values and critical exponents in the 
SYM regime. Fixed points corresponding to a second-order 
phase transition of the system exist in this regime only for 
iVL = 1, 2. 



iVL 


hi 


m« 


A^ 


* 


vl 


Vr 


V 




1 


4.496 


0.326 


5.099 


0.716 


0.142 


0.142 


1.132 


0.786 


2 


3.364 


0.104 


3.643 


0.512 


0.162 


0.325 


1.100 


0.809 



We emphasize that a corresponding scalar 0(27Vl) 
model does not exhibit a fixed-point potential in the SYM 
regime but only in the SSB regime. We conclude that the 



Of course, the full truncation error introduced by the 
derivative expansion is much harder to determine and 
depends on the specific quantity. By analogy with the 
bosonic 0{N) models, we expect the leading critical ex- 
ponent to be our most accurate quantity. On the same 
truncation level, the critical exponent v in 0{N) mod- 
els agrees with the best known value already on the 3% 
level, see [l^. The subleading exponents as well as the 
anomalous dimensions usually are less well approximated 
to this order of the derivative expansion and require more 
refined techniques for a better resolution of the momen- 
tum dependence, e.g., as suggested in p^. [2^. 
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B. The regime of spontaneous symmetry breaking 
(SSB) 

For increasing left-handed fermion number, the scalar 
fixed-point potential must eventually lie in the regime of 
spontaneous symmetry breaking. This is already obvious 
from the structure of the flow equations: for large TVl, the 
Goldstonc-likc fluctuation modes dominate the flow of 
the effective potential (|47p (the A^L-dependent fcrmionic 
contribution in Eq. (j47p is field independent and can be 
dropped). Therefore, the potential flow of our model 
approaches that of an 0(2iVL) model in the limit A^l 
oo. The latter is known to exhibit a Wilson-Fisher fixed- 
point potential in the SSB regime with a nonzero k* > 0. 

We observe the transition of the fixed-point poten- 
tial from the symmetric to the SSB regime already near 
A^L = 3. The properties of the Wilson-Fisher fixed point 
in the analogous bosonic model are, of course, quanti- 
tatively modified by the presence of fermionic fluctua- 
tions, but its basic characteristics are otherwise left in- 
tact. Based on these simple observations, we continue 
with an analysis of the fixed-point structure in the SSB 
regime in the remainder of this section. For the SSB 
phase where the effective potential u is minimal at a 
nonzero value k > 0, we use the expansion (dropping 
the subscript k again for simplicity) 



An / - 



A 



ip-^r = '4(p-'^r+'-^ip-^) 



3 ,~ 



2! 



3! 



(64) 



For the flow of k, we use the fact that the first derivative 
of u vanishes at the minimum, u'{k) = 0. This implies 



= dtu'in) 

^ dtK = — 



dtii'{p)\p^. + idtK)u"{n) 



(65) 



The flow equations in the SSB regime are more involved 
due to additional loop contributions which arise from the 
coupling to the vev. This higher degree of nonlinearity 
inhibits a simple analytical study of the flxed-point struc- 
ture at NLO in the derivative expansion. Instead, we use 
an iterative method, starting at the Wilson-Fisher fixed 
point for the analogous 0(2iVL) model. This fixed point 
can be obtained in our system by setting the Yukawa 
coupling to zero, = 0.^ Starting at this Wilson-Fisher 
fixed point of the reduced scalar system, we can obtain 
a fixed point of the full chiral Yukawa system by nu- 
merical iteration. This confirms that the presence of the 
fcrmions gcnerically shifts the scalar fixed-point values 



^ This offers also the possibihty of a cross check: a comparison 
of our 0(2AfL)-model results to the results which are given in 
Il9ll reveals a very satisfactory precision. The remaining minor 
deviations can fully be attributed to the fact that wc use a simple 
NLO-derivative expansion truncated at As. 



only slightly. However, our chiral Yukawa system rep- 
resents a different universality class, and so the critical 
exponents and the anomalous dimensions are special to 
our system. Numerical results are displayed in Figs.[3][5l 
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FIG. 3: Nonuniversal fixed point values in the SSB regime. 
Left panel: the potential minimum k* in our U(A^l)l®U(1)r 
model (red dots) are compared with those of the 0{2Nh) 
model (blue circles). The negative k* values for A'l < 3 are 
particular to our model and correspond to the fact that the 
fixed-point potential is in the symmetric regime. This transi- 
tion near A^l ^ 3 also causes a kink in the A^l dependence of 
h*^ (right panel). 

In Fig. [3l the fixed-point values of k and are plot- 
ted for different TVl. Note that k* docs not satisfy the 
constraint k > for Ni^ < 3; these negative values are 
therefore unphysical. Of course, this negative branch cor- 
responds to the fixed-point scenario in the symmetric 
regime discussed in the previous section. As expected, 
the fixed-point values of our U(A^l)l®U(1)r model (red 
dots) change only slightly in comparison with the anal- 
ogous 0(2iVL) model (blue circles). Both models were 
investigated with anomalous dimensions and expanded 
up to order ~ in the effective potential. It should 
be stressed that the fixed-point values themselves are 
nonuniversal quantities depending on the details of the 
regulator. 

By contrast, the fixed-point values of the anoma- 
lous dimensions as plotted in Fig. 2] are universal (even 
though slight regulator dependencies can be induced 
by the truncation). The left panel shows 77^ for our 
U(AfL)L'8'U(l)R model (red dots) in comparison with the 
analogous 0(2A'l) model (blue circles). Whereas ry^ in 
both models approaches a common value for large TVl, 
we observe larger differences for smaller TVl which can 
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directly be attributed to the fermionic loop contribu- 
tions. The fermion anomalous dimensions rj^ and rj^ are 
shown in the right panel (purple diamonds and green 
triangles, respectively). For A^l = 1, both anomalous di- 
mensions agree as this corresponds to the left / right sym- 
metric point. For larger A^l, ?7r becomes significantly 
larger, as the massless fermion and boson degrees of free- 
dom contribute to the loop diagrams with a weight ~ Ni^. 
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FIG. 4: Left; in the U(AfL)L®U(l)R model (red dots) in 
comparison with that of the analogous 0(2A^l) model (blue 
circles). Right: (purple diamonds) and-q^ (green triangles) 
in the chiral Yukawa model. 

The critical exponents v and to are shown in Fig. [5l 
Again, the red dots represent the values for our chiral 
U(iVL)L<X)U(l)R model, which is compared with those 
of the analogous 0(2iVL) model (blue circles). The re- 
sults of both models approach each other for large A^l 
as expected but show sizable deviations at smaller TVl. 
We observe a rapid change of the critical exponents near 
TVl — 3, where the effective fixed point potential changes 
from the symmetric to the SSB regime. This strong vari- 
ation is rather natural, as the structure of the fixed-point 
equations varies significantly across this transition. 

We believe that especially the results for the univer- 
sal fixed-point anomalous dimensions and critical expo- 
nents as summarized in Tab. IIIII can provide for first 
benchmarks for these new universality classes. For any 
other nonperturbative tool for chiral fermions near the 
symmetry-breaking phase transition, these universality 
classes can serve as a useful testing ground. As in the 
symmetric regime discussed above, we consider our re- 
sults for the exponent v to be reliable on the few-percent 
level, whereas the subleading exponent as well as the 
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FIG. 5: Critical exponents as a function of A^^l in the chiral 
U(7Vl)l®U(1)r model (red dots) in comparison with those of 
the analogous 0(2A'^l) model (blue circles). The rapid change 
of the exponents occurs near A^l — 3, where the fixed-point 
potential shows a transition from the symmetric to the SSB 
regime. 



TABLE III: Fixed point values and critical exponents in the 
SSB regime. 



AS 



3 
4 
5 
10 
20 
50 
100 



2.718 
2.713 
2.519 
1.452 
0.739 
0.296 
0.148 



0.009 
0.042 
0.079 
0.256 
0.597 
1.612 
3.301 



2.967 
2.954 
2.717 
1.506 
0.752 
0.298 
0.149 



Vr 



0.371 
0.279 
0.204 
0.075 
0.032 
0.012 
0.006 



0.154 
0.125 
0.100 
0.046 
0.022 
0.009 
0.004 



0.487 
0.637 
0.746 
0.913 
0.963 
0.986 
0.993 



0.883 
1.043 
1.124 
1.092 
1.043 
1.017 
1.008 



0.675 
0.678 
0.715 
0.872 
0.942 
0.978 
0.989 



anomalous dimensions might be somewhat less accurate. 



V. CONCLUSIONS 

We have investigated the critical behavior of three- 
dimensional relativistic fermion models with a U(A'l)l x 
U(1)r, chiral symmetry. We have designed this class of 
models to exhibit similarities to the Higgs- Yukawa sector 
of the standard model of particle physics. As a funda- 
ment of this model building, we have classified all pos- 
sible four-fermion interaction terms invariant under this 
chiral symmetry, and also determined the corresponding 
discrete symmetries. We have identified a scalar parity- 
conserving channel similar to the standard-model Higgs 
scalar. For sufficiently strong correlations in this channel, 
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a second-order phase transition into the chiral-symmetry 
broken phase can occur which is a 3d analog of the elec- 
troweak phase transition. Using the functional RG, we 
have computed the critical behavior of this phase transi- 
tion in terms of the critical correlation-length exponent 

the subleading exponent uj, and the fermion and scalar 
anomalous dimensions as a function of Ni^ > 1. 

Whereas the standard model is defined with a fun- 
damental Higgs scalar, we have here started with only 
fundamental fermion degrees of freedom. The result- 
ing Higgs field arises as a scalar bi-fermionic composite 
upon strong fermionic correlations. Therefore, the rela- 
tion between our original fermion model and the result- 
ing Yukawa system is similar to that between top-quark 
condensation models [l^l and the standard-model Higgs 
sector. From a more general RG perspective, however, 
the purely fermionic models can anyway be viewed as 
just a special case of the more general Yukawa models 
j26l . [27! supplemented with nonunivcrsal compositcness 
conditions. As long as the compositeness scale in the 
deep UV remains unresolved, there is no real difference 
between the purely fermionic or the Yukawa-model lan- 
guage. 

For our quantitative results, we have used a consis- 
tent and systematic expansion scheme of the effective ac- 
tion in terms of a nonperturbative derivative expansion. 
Whereas there is an extensive body of circumstantial ev- 
idence in the literature that this expansion is suitably 
adjusted to the relevant degrees of freedom of Yukawa 
systems, a practical test for convergence is problematic 
in the present case. This is because the leading-order of 
the expansion (defined by setting all t^^.l.r = 0) does 
not support the desired fixed point. The latter becomes 
visible only from next-to-leading order on, due to the 
structure of the Yukawa flow. Hence, a straightforward 
convergence test in principle requires a NNLO calcula- 
tion. Instead, the reliability of the results can also be 
verified indirectly: first of all, the derivative expansion is 
based on the implicit assumption that momentum depen- 
dencies of operators do not grow large. This includes the 
kinetic terms, such that self-consistency of the derivative 
expansion requires that the anomalous dimensions satisfy 
^0,L,R ^ 1, as is the case in our calculations. Second, our 
models can always be compared to the purely bosonic 
limit of scalar O(A^) models where the quantitative reli- 
ability of the derivative expansion has been verified to a 
high level of significance. Third, we observe very good 
convergence properties of the polynomial expansion of 
the effective potential which is again a strong signature 
of self-consistency. 

From our classification of all fermionic interaction 
terms compatible with the required symmetries, it is clear 
that the Higgs-like condensation channel is not the only 
possible channel. Aside from vector-like channels, there 
are two further pseudo-scalar channels, cf. Eq. (|32p . and 
further scalar and pseudo-scalar channels in the flavor- 
singlet Fierz transforms of Eqs. (|3ip - ([M|) . In fact, the 
present analysis is a restricted study of a particular con- 
densation process. We expect that the phase diagram of 



the general model is much more involved and might ex- 
hibit a variety of possible phases and corresponding tran- 
sitions. This phase diagram is parameterized by up to 6 
independent couplings being associated with the linearly 
independent fermionic interactions. The calculation of 
the true condensation channel for a given set of initial 
couplings remains a challenging problem. As such a prob- 
lem of competing order parameters is well known also in 
other systems, e.g., in the Hubbard model [1^, H^, [30| . 
the present system can serve as a rich and controllable 
model system. 

A special feature of our model arises in the sym- 
metric regime: here, a fixed point within our trunca- 
tion implies a sum rule for the anomalous dimensions, 
J?!/) + + = 1 = 4 — d. This sum rule is relevant, since 
the underlying balancing between anomalous dimensions 
and dimensional power-counting scaling can be a decisive 
feature of many other models as well. Most prominently, 
the asymptotic-safety scenario in quantum gravity 31 1 
as well as in extra-dimensional Yang-Mills theories [32 1 
requires similar sum rules to be satisfied. In contrast to 
these latter models, the present models for A^l = 1,2 
can serve as a much simpler example for a test of this 
sum rule at a fixed point. A verification of this sum rule 
also by other nonperturbative tools can shed light on this 
important mechanism to generate RG fixed points. 

This is another reason why we believe that the new uni- 
versality classes defined by our models can serve as pro- 
totypes for studies of strongly correlated chiral fermions 
in general and of nonperturbative features of standard- 
model-like chiral symmetry breaking in particular. 
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APPENDIX A: FLOW EQUATIONS AND 
ANOMALOUS DIMENSIONS 



The derivation of the flow of the dimensionless effective 
potential dfU^ (|47p and of the anomalous dimensions 77L 
and r/R (j50|) are explained in the appendix of Con- 
sequently the flow of the coupling constants Xi and the 
flow of the squared mass in the symmetric regime 
and the flow of the dimensionless squared vacuum ex- 
gctation value k in the SSB regime are the same as in 

For deriving the flow of the squared Yukawa coupling 
constant we split the bosonic field into its vev v and 
the deviation from the vev (the relation between the di- 
mcnsionful vev v and the dimensionless squared vev k is 
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V2 



1 

71 



5{p) 



I A</>Kp) + iA^^(p) \ 
A</>?(p) + iA02(p) 



V2 



VA0f-(p) + iA^^^-(p)/ 



(Al) 



We are mainly interested in the Yukawa coupling between the fermions and the Goldstone boson. Thus we use the 
A02 part for the projection. Using the truncation ([5^ and comparing with the Wetterich equation we get 



V25 



't'l-k 



2 #i(p)^A0i(p') 



STr 



dtHvf + Rk] 



p'—p—q—O 



(A2) 



(2.) 

Next, we split {r,:' + Rk) into a propagator part V, which contains only the vev, and a fluctuation part JF, which 



contains the fluctuating fields. Inserting this into Eq. (jA2p the expansion of the logarithm reads 

2 -, / ^\ 3 



In r 



(2) 



Ri 



In 



V 1 



J" 
V 



\n{V) 



T 1 (T 



- - 



(A3) 



Only the term to third power survives the projection. Performing the matrix calculations and taking the supertrace, 
we get 



dthl 



The potential on the right-hand side is evaluated at the minimum ^w^. Using the threshold function and switching 



over to dimcnsionlcss quantities, we end up with Eq. (j49p . 

For the derivation of the flow of rj^, we use the decomposition of (|Aip . Again we use Ai/ij for the projection and 
expand the logarithm as in (jA3p . This time only the quadratic term survives the projection, 



dfZ. 



1 d 



4 V (5A0^(p) 6/\4>\{q) 



STr 



After performing the matrix calculations and taking the mensionless quantities. By use of the threshold functions 



supertrace, we use r/^ 



and switch over to di- (ED), we obtain Eqs. (HOI) in the main text. 
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